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Abstract. We analyze a model of active Brownian particles with non-linear fric- 
tion and velocity coupling in one spatial dimension. The model exhibits two modes 
of motion observed in biological swarms: A disordered phase with vanishing mean 
velocity and an ordered phase with finite mean velocity. Starting from the mi- 
croscopic Langevin equations, we derive mean-field equations of the collective 
dynamics. We identify the fixed points of the mean-field equations corresponding 
to the two modes and analyze their stability with respect to the model param- 
eters. Finally, we compare our analytical findings with numerical simulations of 
the microscopic model. 

1 Introduction 

The concept of active Brownian particles was introduced more than a decade ago 1 1I2I3| and 
was applied to different problems such as for example structure formation in excitable and 
chcmotactic systems [415] or diffusion and swarming in biological systems [617] . A general prop- 
erty of an active system is the emergence of complex behavior in the absence of external forces, 
which may be associated with internal degrees of freedom. On the one hand, there are many 
biological systems which can be interpreted as active particles ranging from individual animals 
[5] to molecular motors [9] . On the other hand, recent reports of artificial systems far from equi- 
librium show characteristics of active Brownian motion |10I11I12"] and demonstrate the general 
applicability of the concept to a wide range of natural phenomena. 

In this work, we will analyze the collective motion (swarming) of active Brownian particles 
with a velocity alignment interaction. Collective motion of living organisms is a fascinating self- 
organization phenomenon which has attracted scientists from various disciplines (13|14ll5] and 
has led also to numerous publication in the field of statistical physics, nonlinear dynamics and 
pattern formation [16I17I18I19I20I21J . Whereas most studies concentrates on simple models of 
collective motion of self-propelled particles with constant speed, only few publications consider 
collective motion of (active Brownian) particles with non-trivial speed dynamics [22 23 24 7 25 . 
In this work we derive in a systematic way mean-field equations for collective motion of active 
Brownian particles. We discuss the analytic results and compare the prediction of the mean 
field theory with numerical simulations. 



2 Microscopic Model 

We consider a one-dimensional system of N active Brownian particles with mass m = 1. The 
evolution of the particle positions Xi and velocities vi is described by the following set of 
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(stochastic) differential equations (i = 1 . . . N) : 

Xi = Vi (1) 

miii — (a - fivf)vi + n(u Sji - v 4 ) + \/2~D& (2) 

The first term on the right hand side is the non-linear friction force, the so-called Rayleigh- 
Helmholtz friction, of active Brownian motion, studied in [26 27 28 . The second term describes 
a velocity alignment interaction with 

1 Ne 

6 3=1 

i.e. with the mean velocity of particles within a finite distance \xj — x,| < e around the focal 
particle i (N e - number of neighbors within e); fj, denotes the alignment strength, which is 
positive fj, > 0. Finally, the last term on the right hand side of Eq. (0) is a white Gaussian noise 
with intensity D and (&) = 0, (&(*)&(*')) = S i:j S(t - £'). 

For a large number of neighbors N e ^> 1 moving with random velocities (disordered state) , 
the mean velocity of neighboring particles vanishes u E: i = 0. In this situation, the mean squared 
stationary velocity of a single particle is given as Vq = (a — (J,)//3. For perfectly ordered motion 
where all particles move with the same velocity (u Sl i = Vi, D = 0), the alignment force vanishes 
and we obtain Vq — a/ '/3. 



3 Mean Field Theory 

We derive a mean field transport theory of a gas of active Brownian particles based on the 
formulation of moment equations for the particle probability distribution. In general, for a sys- 
tem far from equilibrium the probability distribution is not Gaussian, and a correct description 
requires infinitely many moments (see for example |29l30p . Here we perform a closure of our 
moment equations by neglecting temperature fluctuations 9 (Eqs. (|9dp . (jTU)) ). The approxi- 
mation of a non-Gaussian probability distribution by a finite number of moments may lead 
to unphysical behavior, such as negative values or artificial oscillations of the (approximated) 
probability distribution. Therefore, we will later compare the analytical results to numerical 
simulations of the microscopic system. 

The n-th moments of the velocity (v n ) is defined as 

MJx,t) = (v n ) = - [ v n P(x,v 7 t)dv, n>0, (4) 
P J 

where P(x,v,t) is the probability density function describing the probability to find a particle 
at time t, at position x moving with velocity v. The normalization p is the zeroth moment 
which is equivalent to the marginal density 

M (M) = p{x,t) = J P(x,v,t)dv. (5) 

Multiplying the n-th moment with the density and taking the derivative with respect to time, 
we obtain the dynamics of the moments of velocity 



|(M M n) = /^. 



(6) 



We start with an effective single particle description and omit for simplicity the individual 
particle index i. The Fokker-Planck-Equation for a single particle in the mean velocity field u e 
of others reads 

BP 8 8 rP- 
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Inserting Eq. Q in Eq. ((B]) and using \im v ^± oa P(x, v, t) = 0, the terms with partial derivatives 
with respect to v can be partially integrated, yielding 

^(M M„) =-^P (v n+1 ) +np[a (v n ) - (3 (v n+2 ) + fi(u e (v n ^ 1 ) - (v n ))] 

+ n(n- 1) Dp (v n ~ 2 ). (8) 

We rewrite velocity of the focal particle as a sum of the local velocity field u(x,t) plus some 
deviation Sv: v = u + Sv. Furthermore, we assume (Sv ) = for odd exponents I (I = 1, 3, 5, . . . ). 
Thus we obtain for the moments (up to I — 4): 

(v) = u , (9a) 

(v 2 ) =u 2 + T , (9b) 

(v 3 ) = u 3 + 3uT , (9c) 

(v 4 ) = m 4 + 6 u 2 T + T 2 + 9 . (9d) 



Here, T is the mean squared velocity deviation T = (Sv 2 ), which we will refer to as the tem- 
perature of the active particle gas, whereas 9 is the average of the mean squared temperature 
fluctuations defined as 

9 = (((v-u) 2 ~Ty) = (5v 4 )-T 2 . (10) 

Now we can insert the Eqs. © in Eq. ([5]). Considering the dynamics up to n — 2, after some 
calculus we arrive at a set of three coupled partial differential equations for the evolution of the 
density p(x,v,t), the mean velocity field u(x,v,t), and the temperature field T(x,v,t): 

^P = "^(P«) (Ha) 
du ^ 9 fl ( *±vr\ j- ( , dT T dp 

_ + u _ u = au _p u ( u + 3r) +M ( Ue - u) -__-_ (nb) 



Kf + ^)= {a -» ) T- P TW + T ) -^ + D-T d £ 



(11c) 



Let us consider for simplicity an isotropic system with vanishing gradients in mean velocity u 
and temperature T, which is a good approximation of our systems at high particle densities. In 
this case, the local velocity in the velocity alignment force equals the constant mean velocity 
field across the system u s = u, and we end up with the following two ordinary differential 
equations for the temporal evolution of u and T: 

^au-/3ii (u 2 + 3T) (12a) 
=(a - n) T - T (3u 2 +T)-(39 + D. (12b) 

In order to obtain a closed system of equations, we neglect the temperature fluctuations by 
setting 9 = in Eq. (|10p . which is a reasonable assumptions at small noise intensities. Thus, 
the above differential equations constitute a two-dimensional dynamical system, with 6 fixed 
points (stationary solutions) in the (u, T) phase space, which can be analyzed by means of 
linear stability analysis. 
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Fig. 1. Bifurcation diagram of the mean velocity u vs noise intensity D as predicted from the mean 
field theory for different velocity alignment strengths (a = ft = 1) : (a) /i < 2a/3 (/i = O.f), (b) 
2a/3 < n < 10a/3 (fi = 0.7) and (c) n > 10a/3 (fi = 5.0). (d) Phase diagram with respect to velocity 
alignment \x and noise intensity D. 



The stationary solutions {du/dt = dT/dt = 0) for u and T read: 

u h2 = 0, (13a) 

Tl ' 2 = 2? (13b) 



10a - 3 (ji - ^{2a + ii) 2 - 32 (3D) 
u 3A = ± J (13c) 

T 3 =T 4 = — (13d) 

J 10a -3U+ v /(2a + n) 2 - 32/3D) 

u " = 4^ (13c) 

n,-n,- *-^vyy- i °ff p . m 

The kinetic temperature has to be positive, therefore, for u — 0, T\ (positive square root) is 
the only physically reasonable solution. 

The first solution with vanishing mean velocity u — describes a disordered phase. For 
D = 0, the temperature T = ^j^- = Vq equals the square of the stationary velocity of individual 
particles. The kinetic energy of all particles completely consists of fluctuations, no systematic 
translational motion occurs. 

The second pair of solutions corresponds to translational modes which are stable below a 
critical noise intensity. The two solutions correspond to translational motion with positive or 
negative velocity u, thus, to a collective motion of the particles to the left or right. Without 
noise, — 0, and the stationary mean velocity reduces to 113,4 = ±^/a~/P- Increasing the 
noise rises the kinetic temperature, and results in a decrease of the mean speed \u\. 
The last solution pair describes unstable modes, for which with increasing noise intensity D 
the temperature decreases and the mean speed increases. 

For low velocity alignment strength fi < 2a/3, the disordered phase is always a stable 
solution; for fi > 2a/3, the linear stability analysis of the mean-field equations predicts the 
existence of a critical noise intensity 

- (14) 

which determines the stability of the disordered solution. Starting from large noise intensities 
where the disordered solution is stable and decreasing the noise below £>i iC rit, we observe a 
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Fig. 2. (Color online) Comparison of the stationary solution obtained from simulation with theoretical 
prediction from the mean field theory for different velocity alignment strengths: fj, — 0.4 (a,d), /j, — 0.67 
(b,e) and /i = 5.0 (c,f) and different friction coefficients: /? = 0.25 (a,b,c) and ft = f .0 (d,e,f). Other 
parameters used: particle number TV = 8192, simulation domain L = 500, velocity alignment range 
e — 50 and a — 1.0. The initial conditions were either the disordered state (black filled squares) or the 
ordered state (red circles). Solid (dashed) lines show the stable (unstable) stationary solutions of the 
mean-field equations. 



pitchfork-bifurcation, and the disordered phase becomes unstable. Depending on the value of 
/j,, the pitchfork-bifurcation is either sub- or super-critical. For fi < 10a/3, the disordered 
solution becomes unstable through a collision with the two unstable translational solutions, 
whereas for fi > 10a/3 no unstable translational solutions exist and the disordered solution 
becomes unstable directly through the appearance of the two stable translational solutions (see 
Fig- HI- Thus, for \i > 2a/3 and D < £>i iCr it, only the translational solutions 1*3.4 are stable. 

For /i < 10cv/3 there exists a second critical noise intensity which determines the stability 
of the ordered phase (translational solutions, u^O). Above the critical noise intensity 

D 2crit = ^i# (15) 

■ 32/3 v ; 

all translational solutions become unstable through a saddle-node bifurcation (Fig. Q] a,b) . 

In order to test our analytical results, we performed numerical simulations of the microscopic 
model with periodic boundary condition. Due to the symmetry of the translational solutions 
«3 = — U4, we distinguish the disordered phase and the ordered (translational) phase by mea- 
suring the global mean speed in our simulations: 



<M> = 



1 N 



(16) 



Here, (•) denotes temporal average after the system has reached a stationary state. In order 
to analyze the stability of the (dis)ordered phase, the simulations were performed with two 
different initial states: perfectly ordered with u(t = 0) = u^(D = 0) and perfectly disordered 
state with u(t = 0) = 0. Each simulation was run until the system reached a stationary state 
but at least for t = 2000 time units with a numerical time step At = 0.01. 

The stationary speed of the ordered phase vs noise intensity obtained from numerical simu- 
lations with ordered state initial condition are in a good agreement with the theoretical predic- 
tions from the mean field theory. Whereas simulations with disordered initial condition reveal 
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Fig. 3. (a-d) Simulation snapshots showing the breakdown of the disordered solution. Each point 
corresponds to a single particle. Starting from a perfectly disordered initial state ((a), t — 10), we 
observe a local symmetry breaking at x ~ 400 ((b), t = 160) leading to a local velocity alignment of 
particles. The region of alignment spreads through the system ((c), t — 210) until a stationary ordered 
state is reached ((d), t — 350); Simulation parameters: a = 1.0, /3 = 1.0, fi = 0.4, D = 0.01, N = 4096, 
L = 500, e — 10. (e) Stationary temperature for non-interacting particles (/i = 0) vs noise intensity D 
obtained from the mean field theory (solid line), from direct calculation of the second moment Eq. (|17|l 
(dashed line) and from numerical simulations (symbols) for a = /? = 1. 



an unexpected instability of the disordered solution. At finite \x the numerical simulations show 
that at intermediate D the disordered solution u = becomes unstable via a spontaneous 
symmetry breaking as shown in Figs. [2] and [3^,-d, which is not predicted by the mean field 
theory. 



4 Discussion and Summary 



We confirmed the instability of the disordered solution in numerical simulations for interme- 
diate noise strengths for different particle numbers (N — 4096, 8192, 16384) and obtained the 
same result for all values of N. This suggests that it cannot be simply dismissed as a pure fi- 
nite size effect. Other possible sources for this discrepancy might be the restriction to spatially 
homogeneous solutions of the mean-field equations, which neglects density fluctuations, or the 
closure of the moment equation hierarchy by neglecting temperature fluctuations 9. The latter 
is consistent with the observation that the deviation does not appear at low D, where tempera- 
ture fluctuations are very small. In fact, for non-interacting particles with Rayleigh-Helmholtz 
friction we can calculate directly the second moment from the velocity distribution |26j . For 
P = 1 it reads 




with /„ and K n being the modified Bessel functions of the first and second kind respectively. 
The result of Eq. (| 1 7[) corresponds directly to the temperature T for the disordered state in 
the limit /i = 0. Due to the non-linearity of the friction function, the temperature does not 
increase monotonically with D as predicted by the mean field theory but exhibits a minimum 
at intermediate noise intensities as shown in Fig. [3^. It can be seen from Eq. (TT2"]) that the 
explicit consideration of finite temperature fluctuations (9 > 0) leads to a decrease in T, which 
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is consistent with Eq. (IT71) . This in turn decreases the stability of the disordered state. The 
extension of the mean field theory to higher orders would account for this effect at the expense 
of the analytical tractability of the mean field solutions. 

In summary, we have formulated a one-dimensional model of active Brownian particles with 
non-linear Rayleigh-Helmholtz friction interacting through a local velocity alignment. Starting 
from the microscopic Langevin description, we derived a set of mean-field equations via the 
corresponding Fokker-Planck Equation. By neglecting temperature fluctuations, we obtained 
a set of coupled partial differential equations for the density, velocity and temperature fields. 
For simplicity, we restricted the analysis here to the spatially homogeneous situation. We have 
identified mean-field solutions corresponding to an ordered mode (collective motion) and to the 
disordered mode with random particle velocities and analyzed their stability. A comparison of 
our analytical results with numerical simulations of the microscopic model confirms our findings 
for the ordered phase, but shows deviations with respect to the stability of the disordered 
solutions. 

P. R. would like to thank M. Kostur for the introduction to numerical simulations on 
Graphical Processing Units using CUDA which allow a considerable speed up of numerical 
simulations [3T]. U. E. would like to thank B. Eckehardt for fruitfull discussions on active 
Brownian particles. Last but not least, we would like to express our gratitude to L. Schimansky- 
Geier for his scientific input in our research and the time we went together; we wish him all 
best in the future. 
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